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Abstract —This study investigates wireless information and 
energy transfer for dual-hop amplify-and-forward full-duplex 
relaying systems. By forming energy efficiency (EE) maximization 
problem into a concave fractional program of transmission 
power, three relay control schemes are separately designed to 
enable energy harvesting and full-duplex information relaying. 
With Rician fading modeled residual self-interference channel, 
analytical expressions of outage probability and ergodic capacity 
are presented for the maximum relay, signal-to-interference-plus- 
noise-ratio (SINR) relay, and target relay. It has shown that 
EE maximization problem of the maximum relay is concave for 
time switching factor, so that bisection method has been applied 
to obtain the optimized value. By incorporating instantaneous 
channel information, the SINR relay with collateral time switch¬ 
ing factor achieves an improved EE over the maximum relay in 
delay-limited and delay-tolerant transmissions. Without requiring 
channel information for the second-hop, the target relay ensures a 
competitive performance for outage probability, ergodic capacity, 
and EE. Comparing to the direct source-destination transmission, 
numerical results show that the proposed relaying scheme is 
beneficial in achieving a comparable EE for low-rate delay-limited 
transmission. 

Index Terms —Energy harvesting, wireless power transfer, 
amplify-and-forward relay, full-duplex relay, relay gain control. 


1. Introduction 

Energy harvesting (EH) has emerged as a promising ap¬ 
proach to prolong the lifetime of energy constrained wire¬ 
less communications (Hl-El- Through harvesting energy from 
natural sources (e.g., solar, wind, thermoelectric effects or 
other physical phenomena), periodic battery replacement or 
recharging can be alleviated in green communications Q-lEl. 
However, EH from natural sources is not typically applicable 
for devices with small sizes or located in physically limited 
places in, in. With the capability to harvest energy from 
ambient radio-frequency (RE) signals, simultaneous wireless 
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information and energy transfer (SWIET), also known as si¬ 
multaneous wireless information and power transfer (SWIFT), 
provides a more promising way for green wireless commu¬ 
nications to function in environments with physical or other 
limitations a-im. 

The pioneering work on SWIET can be traced back to 
16) and 0, where the fundamental tradeoff between capacity 
and energy was studied for point-to-point communications. 
Eollowing the assumption that an ideal receiver is capable of 
observing information and extracting energy from the same 
received signal, SWIET has been extended to multi-antenna 
systems ca, ina, multiuser systems M, and bi-directional 
communication systems US- However, as discussed in il, a 
practical circuit for EH from the RE signal can hardly decode 
the carried information from the same signal. Therefore, two 
practical receiver architectures, namely, time switching (TS) 
and power splitting (PS), are proposed in They are now 
widely adopted in various wireless systems, such as multiple- 
input multiple-output (MIMO) systems iTEl . orthogonal fre¬ 
quency division multiplexing systems |[T7]| . and cellular sys¬ 
tems ll9lFI . 

In parallel with the aforementioned studies that mainly deal 
with single-hop scenarios, employment of cooperative relays 
to facilitate RE EH and information transfer in energy-efficient 
cooperative or green networks has also drawn significant 
attention lH, ifTOl . ifTTIl . |[T^ . Relay-based SWIET not only 
enables wireless communications over long distances or across 
barriers, but also keeps the energy-constrained relays active 
through RE EH. The authors of ifTTIl designed and analyzed 
the TS and PS relaying protocols for amplify-and-forward (AE) 
relaying systems, and then extended the results to an adaptive 
TS relaying protocol llT9l . The throughput of the TS and 
PS relaying protocols for decode-and-forward (DE) relaying 
was investigated in Gol. Several power allocation schemes 
for EH relay systems with multiple source-destination pairs 
were studied in ||20l. The outage and diversity of SWIET 
for cooperative networks with spatially random relays were 
investigated in EH and the distributed PS-based SWIET was 
designed for interference relay systems ll22l . More recently, 
antenna switching and antenna selection for the SWIET re¬ 
laying systems have been investigated in ll^ and ll24l . re¬ 
spectively. Nevertheless, all these studies are limited to half- 

^Note that antenna switching can also be applied in the multiple-antenna 
case (13. 
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duplex relaying (HDR) mode. Since the source-to-relay and 
relay-to-destination channels are kept orthogonal by either 
frequency division or time division multiplexing, significant 
loss of spectral efficiency occurs in the HDR mode. As an 
alternative, full-duplex relaying (FDR) has drawn considerable 
attention Ga-iEol. Since FDR requires only one channel for 
the end-to-end transmission, a significant improvement in the 
spectral efficiency over the HDR can be achieved. 

So far, a few studies have been conducted for SWIFT in 
FDR systems. In ll27]| . the throughput of the TS relaying 
protocol has been analyzed for SWIFT FDR systems, in 
which the EH relay is operated cooperatively. In practice, 
since the relay node suffers severe self-interference from its 
own transmit signal, FDR operation is difficult to implement. 
For example, self-interference of more than 106 dB has to 
be suppressed by a femto-cell FDR base-station to achieve 
the link signal-to-noise-ratio (SNR) equals to that of an HDR 
counterpart ED. For systems that require a higher transmis¬ 
sion power, more self-interference suppression is needed ll^ . 
Therefore, MIMO has been employed at EH FDR node to 
suppress self-interference ll3^ . By employing EH relay in 
the second time phase for the conventional two-phase AF 
HDR systems, a self-interference immunizing full-duplex relay 
node was proposed in ES, which can transmit information 
and extract energy simultaneously via separated transmit and 
receive antennas. Note that all the above researches of EH 
FDR are conducted to maximize throughput, whereas energy 
efficiency (EE) has not been investigated. Another challenging 
problem for SWIET is to determine the EH parameters, i.e., 
the TS factor and PS factor for the TS and PS relaying 
protocols, respectively. Determining of the TS factor affects not 
only the relay-harvested energy but also the effective relaying 
transmission time in the TS relaying protocol. Compared with 
the PS relaying protocol, the TS relaying protocol is more 
practical because of its simplicity. With statistical channel 
state information (CSI), the numerical optimizations of the 
TS factor have been presented in ifTH and lIZTll . In delay- 
limited and delay-tolerant transmissions, the instantaneous CSI 
can also be applied to optimize the EH parameter, as well as 
the instantaneous CSI-aided transmission power control in the 
EDR systems ||28l, 1^ . 

Motivated by these previous studies, we focus on wireless 
information and energy transfer for an energy-constrained dual¬ 
hop EDR system, in which the source node has a reliable 
green power supply, whereas the TS-based AE relay node is 
powered via EH from the source-emitted RE signall. Since 
self-interference deteriorates EDR performance seriously, relay 
gain control with SWIET becomes more critical than ever 
before. Compared with existing works, some distinct features 
of our study are highlighted here. In ||33l, the effective infor¬ 
mation transmission time is the same as that of HDR systems. 

^This setting has a number of potential applications in green wireless 
network, e.g., when the intermediate node is energy-selhsh or lacks an energy- 
supply, or the direct link from the source to destination is blocked by a barrier 
while the relay has to be placed on a site without a fixed power supply. 


Therefore, the spectral efficiency improvement is minimal 
compared with that of the EDR systems. In our study, we focus 
on improving EE by applying full-duplex information relaying, 
i.e., the relay receives and forwards the source information to 
the destination simultaneously. In IZTll . the authors designed 
wireless information and energy transfer for EDR systems 
with an aim to maximize throughput. Eor full-duplex radios, 
experiments have shown that residual self-interference (RSI) 
contains specular component and RSI channel is characterized 
by Rician fading with its iT-factor depending on the amount of 
direct-path suppression El- Similar trends for RSI characters 
can also be observed in El- Since EE is a major concern in 
the design of SWIET EDR, we focus on the impact of this 
experimental verified RSI channel model on EE performance 
by optimizing TS factor, relay gain, and source transmission 
power. 

In this paper, three relay control schemes, namely, the 
maximum relay, signal-to-interference-plus-noise-ratio (SINR) 
relay, and target relay, are investigated for the three trans¬ 
mission schemes of instantaneous transmission, delay-limited 
transmission, and delay-tolerant transmission. The contribu¬ 
tions of this study are summarized as follows: 

• Considering that EDR RSI channel follows Rician fading, 
we present analytical expressions for the throughput for 
the three relay control schemes in different transmissions, 
such as instantaneous transmission, delay-limited trans¬ 
mission, and delay-tolerant transmission. Specifically, an¬ 
alytical expressions for the outage probability are derived 
for these three relay control schemes in delay-limited 
transmission, while analytical expressions for the ergodic 
capacity are derived for the maximum relay and SINR 
relay in delay-tolerant transmission. 

• We show that EE maximization with given TS factor and 
relay gain is a concave fractional program of transmis¬ 
sion power, so that EE monotonically decreases with an 
increasing transmission power in instantaneous and delay- 
tolerant transmissions for practical transmission poweijj. 
In delay-limited transmission, we show that a relatively 
higher outage probability, as well as correspondingly a 
relatively lower transmission power, may achieve a higher 
EE. When the source-destination direct link is not avail¬ 
able, the numerical results also show that SWIET EDR 
is beneficial in achieving EE for low-rate delay-limited 
transmission. By enhancing the performances of outage 
probability, ergodic capacity, effective relaying time, and 
corresponding EE, the SINR relay is shown to be more 
competitive than the maximum relay. Without requiring 
CSI of the second-hop, the target relay also achieves a 
competitive outage probability, ergodic capacity, and EE. 
We also reveal that the worst performances of outage 
probability, ergodic capacity, and EE are obtained by the 
relay placed midway between the source and destination. 

• It has shown that the EE maximization problem of the 

^Transmission power is not very small, e.g., not less than 0 dBm. 
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maximum relay is concave for all possible values of the 
TS factor. Thus, the low-complexity bisection method can 
be applied to obtain the optimized TS factor. The closed- 
form TS factors have been designed for the SINR relay 
and target relay to maximize the e-SINR and achieve 
a target e-SINR, respectively. With the obtained closed- 
form TS factors for the SINR relay and target relay, 
instantaneous CSI has been employed to improve the EE 
in delay-limited and delay-tolerant transmissions. 

The rest of this paper is organized as follows. Section II 
describes the system model of the considered EDR system and 
formulates the EE optimization problem. Section III presents 
the three relay control schemes. The analytical results of the 
throughput are presented in Section IV. Section V presents 
numerical results and discusses the system performances of 
our proposed relay control schemes. Einally, Section VI sum¬ 
marizes the contributions of our study. 

II. System Model 

In this paper, we consider a wireless dual-hop EDR system, 
in which a source node intends to transfer its information to the 
destination node. Due to physical isolations or environmental 
limitations between the source and destination, a cooperative 
relay is employed to assist information transmission from the 
source to the destination. The cooperative relay is assumed 
to be an energy-selfish or energy-constrained device such 
that it needs to harvest energy from the source-emitted RE 
signal to forward the source information to the destination. 
Eor simplicity of implementation, the AE relaying scheme and 
TS transceiver architecture are chosen at the relay node. The 
channel coefficients from the source to the relay and from the 
relay to the destination are denoted by hi and / 12 , respectively. 
hi and /12 are assumed to be frequency non-selective and 
quasi-static block-fading, following a Rayleigh distribution. 
The means of the exponential random variables \hi\‘^ and 
|/i 2 p are denoted by Ai and A 2 , respectively. According to 
the experimental results of (341, the RSI at a full-duplex 
node contain specular components such that the RSI channel 
coefficient can be modeled as a Rician random variable ho. 
Since both experiment and analysis have shown that RSI 
cannot be eliminated completely, this paper considers only the 
case of \ho\‘^ > 0 (341 . (36l . 

The framework of the TS relaying protocol is illustrated 
in Eig. 1(a), in which each time block T is divided into 
two phases. Denoting the TS factor by a (0 < a < 1), 
we use the first phase assigned with a duration of aT for 
energy transfer from the source to the relay. The second phase 
assigned with the remaining duration of (1 — a)T is used for 
full-duplex information relaying via the dual-hop channel. The 
relay-received RE signal in the two time phases are sent to 
the EH receiver and full-duplex transceiver, respectively, as 
illustrated in Eig. 1(b). Since the relay node does not transmit 
during the first time phase, self-interference is not introduced 
during the EH period. The harvested energy at the relay is 


T 


Energy Harvesting 

Full-Duplex Relaying 
Tx/Rx 

◄--► 


aT {\-a)T 


(a) 



(b) 


Eig. 1. (a) Illustration of the key parameters for EH and EDR at the relay, 

(b) Block diagram of the relay transceiver. 


given by 

S = (1) 

where 77 is the energy conversion efficiency depending on the 
rectifier circuit, £1 is a path loss factor, di is the distance 
between the source and relay and it has been normalized with 
respect to the reference distance, m is the path loss exponent, 
and Ps is the source transmission power. In this paper, the 
path loss factor is defined as Ci = AsC, where As is the 
source transmit antenna gain and C is the measured path loss 
at the reference distance. Given the relay-harvested energy, the 
maximum transmission power at the relay is expressed as 

= = ( 2 ) 

where /u = 

In the EDR mode, the relay concurrently receives the signal 
yr{t) and transmits the signal Xr{t) on the same frequency. As 
depicted in Eig. 1(b), the full-duplex transceiver down-converts 
the received RE signal to the baseband, processes the baseband 
signal, and up-converts the processed baseband signal. In Eig. 
1(b), h}a\^) (^) the narrow-band Gaussian noises 

introduced by the receive and transmit antennas, respectively. 
In addition, (t) and (t) are the baseband additive white 
Gaussian noises (AWGNs) caused by down-conversion and 
up-conversion, respectively (H. Eor simplicity, the equivalent 
baseband noise composing both (t) and (t) is modeled 
by the zero mean AWGN ria^ (t) with the variance cr^, and the 
equivalent baseband noise composing both (t) and (t) 
is modeled by the zero mean AWGN ric^ (t) with the variance 
cr^. Therefore, the overall AWGN at the relay node can be 
modeled as the zero mean AWGN nr{k) = ria^k) + ric^ {k) 
with the variance + At the relay node, the sampled 
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baseband signal is given by 

Vrik) = ^J~^^his{k) + ArhoXr{k) + nr{k), (3) 

where k denotes the symbol index, s{k) is the sampled and 
normalized information signal from the source, Xr{k) is the 
sampled signal of Xr{t), and Ar is the relay transmit antenna 
gain. Note that the second term on the right side of (O is the 
RSI. 

Using the harvested energy, the relay amplifies the received 
signal by a relay gain f3. Then, the transmitted signal at the 
relay can be expressed as 

Xrik) = ^/Pyr{k - r), (4) 


where r > 1 is the processing delay at the relay. By recursively 
substituting ^ and dH, we have the following expression for 
the transmitted signal at the relay: 

oo 

Xr{k) = ^^{Arho^)^ 

1=1 

X his{k — It) + nr{k — Ir)^ . (5) 

The sampled received signal at the destination is given by 


yd{k) = ^J^h 2 Xr{k) + nd{k), (6) 

where ^2 is the distance between the relay and destination, 
which is also normalized with respect to the reference distance, 
C 2 = ArC is a path loss factor, and nd{k) is the AWGN with 
the zero mean and the variance aj. By substituting ^ into 
we have 


00 

yd(k) = Y. {ArK^/pAs{k - It) 

1 = 1 
00 

{Arhp^)^ ^Tirik-lr) +nd{k).{l) 


i=i 


In the following, we derive the end-to-end signal power 
under the condition of employing cooperative non-oscillatory 
relays. By assuming that all the signal and noise samples 
are mutually independent, we calculate the relay transmission 
power Pr =K{\xr{k)\‘^} from ® as 




1=1 


0 




( 8 ) 


1-Ar\ho\‘^0 ■ 

To prevent oscillation and guarantee a finite relay transmission 
power, the relay gain is limited by 


0< 


1 


(9) 


The actual relay transmission power should be less than or 
equal to the maximum transmission power, i.e.. 


<Pe- 


( 10 ) 


When Q and ® are substituted into (fTOl) . the relay gain under 
the maximum relay transmission power is limited by 


P< 


_ y _ 

^ TsR + I ^0 P 


( 11 ) 


where the channel SNR of the first-hop link is defined as 
7sr — • ^1 symbol index k, the destination node 

can employ any standard detection procedure to decode the 
desired signal s{k — r), and the rest of the received signal 
components act as interference and noise. Again based on the 
assumption that signal and noise are independent of each other, 
the received signal power at the destination is calculated from 
® as E{\yd{k)\‘^} = + cr^, which 

can be further evaluated as follows, comprising of the desired 
signal power, RSI power, and noise power: 


E{\yd{k)f} = £i£2£r™£2“™/3n|/iinii2| 


desired signal power 


■C2d2"^0\h2\ 


RSI power 


M • 


( 12 ) 


noise power 

Based on (O, the e-SINR at the destination is given by 

7sr7rd 


7srd 


7sr//^ + 7rd + (7sr + 1)7rd l /^- A ^ lhol '^ 


(13) 


where the channel SNR of the second-hop link is defined as 

^ A = CiC2Ps\hi\^\h2\^ 

In this paper, the EEs of instantaneous transmission, delay- 
limited transmission, and delay-tolerant transmission are con¬ 
sidered. The EE is defined as the number of bits successfully 
conveyed to the destination per Joule consumed energy and is 
given by 


%ff{Ps,a,0) = BR^{Ps,a,0)/Ps, (14) 

where B is the system bandwidth and i?j^(Ps, a,/3) is the 
throughput of the transmission scheme M. Let M = I, DL, and 
DT represent instantaneous transmission, delay-limited trans¬ 
mission, and delay-tolerant transmission, respectively. Their 
corresponding throughputs are respectively given by 

R0Ps,a,P) = (1 - Q;)log 2 (l+ 7 ), (15a) 

R^dPs,a,0) = (l-a)(l-Pout)R, (15b) 

R^dPs,a,0) = (l-a)C'E, (15c) 

where Pout = Pr(7 < 7th) is the outage probability, R = 
log 2 (l + 7th) is the fixed transmit rate, = E{log 2 (l + 
7 )} is the ergodic capacity, and 7 th is an e-SINR threshold 
for correct data detection at the destination. The design goal 
is to maximize the EE by optimizing the control parameters 
{Ps , ct,/3}. The optimal control parameters {P*,ct*,y5*} can 




















5 


be obtained by solving the following optimization problem: 


= arginax?7^„(P5,Q;,/3) (16a) 

a,p 

subject to 0 < Ps < Pmax, (16b) 

0 < a < 1, (16c) 

(I6d) 


Since has a very complicated express with respect to 
{Ps,Q(,/ 3} and the information of {a,/3} is unknown, the 
optimal control parameters {P*,Q(*,/3*} can hardly get a 
closed-form solution. By verifying the second derivative of 
7srd with respect to P^, it can be shown is a 

concave function when /3 satisfies the non-oscillation condition 
dg. Since log 2 (l + 'Jsrd) and 1 - Pr( 7 sj,D < 7th) in (15) 
are monotonically increasing with 7gp^D, log 2 (l +7srd) 

1 — Pr(7max < 7srd) ^1^0 concave for all possible 
values of Pg. Thus, the numerator of r]^^^{Ps) is a concave 
function for any {a, f3} satisfying (I16cb and (I16db . Note that 
the denominator of {Ps), Ps, is convex respect to Pg. Then, 
for any {a^/3} satisfying (I16cb and (I16dl) . the optimization 
problem (HSi becomes a concave fractional programing of 
Pg, which can be solved exactly by Dinkelbach’s method 
lEl. The optimization problem can be reduced to optimizing 
{q(,/3}, so that the system design becomes mathematically 
tractable. Since Pout and Ce generally require only statistics 
of CSI, the appearance of the term (1 — ct) in (I15bl) and (I15cl) 
intrinsically implies that a is optimized statistically for the 
delay-limited and delay-tolerant transmissions. However, as we 
will show lately, instantaneous CSI can also be employed to 
optimize a for delay-limited and delay-tolerant transmissions. 
In the following, all the three relay control schemes and 
their EE-related performances are presented for the considered 
transmission schemes, respectively. 

III. Relay Control Scheme 

In this section, we investigate how to compute the relay 
gain and TS factor for the three relay control schemes, namely, 
maximum relay, SINR relay, and target relay. In designing the 
relay control scheme, we assume that perfect knowledge of 
CSI is available. 


A. Maximum Relay 

A simple and popular relay control scheme is involved in 
setting the relay gain at the maximum relay transmission power 
Eg, ED, EBl- In contrast to these works, our study considers 
the maximum relay targeting at maximizing the EE in the 
presence of RSI. Eor a given a in the range (0,1), the relay- 
harvested energy and maximum relay transmission power are 
determined, so that the relay gain is given by the following 
according to (HUl: 

^ 1 + 7 sr+ 7 sr' 4 -I'»oI''’ 


which guarantees that (O holds. Substituting (fTTl) into (fTSl) . 
the e-SINR achieved by the maximum relay is given by 


7sh+(^7sh^.|/»oP + i)(^7r,+i) ■ 

Now, the instantaneous throughput can be evaluated by substi¬ 
tuting ([18]) into (I15al) . 

The appearance of the term in ([13 and (fTSl) 

indicates that the system performance of the maximum relay is 
affected by RSI. Since is a Rician variable, |/ioP follows the 
non-central chi-squared distribution and its probability density 
function (PDE) is given by 




A {K+l)e- 


(R+i)« 




where =E{\ho\‘^} is the average RSI channel gain and K 
is the Rician P-factor. 

Proposition 1: The outage probability achieved by the max¬ 
imum relay is given by 


OO 

1 — ^ J J f\ho\‘^{w)e dzdw(20si) 

w=0z=- 


1 

«!- J fiho\^(w)pKi{p)e~^dw, (20b) 
0 


where a = jCiPgdi^d^a‘^jth{^-\-fJ^Arw), b = 
c = ClC 2 Pgti{l - p^ythArw), d = CiC 2 Psd^a‘^p^th, p = 
V CX1X2 ’ first-order modified Bessel function of 

the second kind ||39j Eq. (8.432)], and the approximation in 
the step (a) is achieved in the region of high SNR values. 

Proof: A proof is provided in Appendix A. ■ 

Proposition 2: The ergodic capacity achieved by the maxi¬ 
mum relay is given by 


CX) 

Y 1 [ 4 ^ Al A2(l + AtAr-'R’) 

'B-hr2 J /|.,|2(^)W^2 

0 


1,0 ) 


dw 


00 ( d'pd^aiKiK+i) y +1 

if In 2 (n!)2 

n=0 

—1—n ^ 

\CiC 2 ri‘^PsXiX 2 Arcrl 0,-l-n,-l-n,-ny ^ 


( 21 ) 


where is the Meijer G-function |[^ Eq. (9.301)]. 

Proof: A proof is provided in Appendix B. ■ 

When the maximum relay is applied, the optimized TS factor 
can be obtained by solving the following optimization problem: 


a* = arg max (a) 
subject to 0 < Of < 1. 
By verifying the sign of we have 


( 22 ) 

> 0 when 


0 < a < 


_ V Tsr + ^_ 

Vo'sr+I+V'IsrT'rd 
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and < 0 when 


_V _ 

+ ^ + '^SR'^RD I^oP 


< a < 1. 


Thus, 7 max is concave for all possible values of a. Since 
log 2 (l +7max) and 1 - Pr(7max < 7th) in (15) are monotoni- 
cally increasing with 7 max, log 2 (l+ 7 max) and 1 -Pr( 7 max < 
7th) are also concave for all possible values of a. Based 
on the fact that the coefficient (1 — a) in (15) does not 
change the convexity of inherited from log 2 (l + 7max) 
and 1 — Pr(7max < 7th), it can be concluded that is 

concave. Then, the optimized ct* can be obtained by solving 

dll (cki) 

the equation —^— = 0. Given the complicated expression 

dll f Cki ^ 

for —= 0 , the closed-form solution is difficult to obtain. 
However, since a unique global maximum of exists, 

the optimized TS factor can be obtained by applying the 
bisection method with the complexity of 0 (log ^), where 5 
is the searching precision mni. Since Pout and involve 
only statistics of CSI, the maximum relay takes no advantage 
of instantaneous CSI in the delay-limited and delay-tolerant 
transmissions, respectively. 


B. SINK Relay 

Since the EE is affected by both the e-SINR and relaying 
transmission time when the source transmission power is 
fixed, TS factor can be optimized to maximize the e-SINR 
and relaying transmission time jointly or separately. In this 
part, a relay control scheme is proposed to maximize the e- 
SINR with the aid of instantaneous CSI. Compared with the 
maximum relay which searches TS factor to optimize the 
e-SINR and relaying transmission time simultaneously, the 
proposed scheme only maximizes the e-SINR and we call it 
the SINR relay throughout the paper. 

According to (O, given the received signal power at the 
destination node as a function of the relay gain, the desired 
signal power is linear, but the RSI power is nonlinear. Conse¬ 
quently, increasing the relay gain can increase the RSI power 
faster than a desired signal power and lead to a reduced e- 
SINR. We can show that (O has a single maximum point for 

G (0, setting the derivative of ([13]) equals 

to zero, the relay gain of the SINR relay that maximizes the 
e-SINR is obtained as 


relay-harvested energy is just enough to implement the relay 
gain /3sinr. Therefore, by solving Psinr = /^max for any a 
(0 < a < 1), we can provide the TS factor for the SINR 
relay as 

Since no redundant energy has been harvested by employing 
asinr, the relaying transmission time (1 —asinr)P is longer than 
those of other as satisfying a > asinr- The TS factor asinr can 
be computed at the destination, or at the relay locally when 
the relay can access the global CSI. 

Proposition 3: The outage probability achieved by the 
SINR relay is given by 


out 


l_j_ f _A±K 

Al j ~i- + K+Ar- 

jO^iPs 


-FC.Ar ctq P 

(26) 


where p = 




7th (l+7th) Jm ‘2 




/: 2 A 2 (/:iPa^-d^crp 7 ^h) 


Proof: A proof is provided in Appendix C. ■ 

Proposition 4: The ergodic capacity achieved by the SINR 
relay is given b>EI 


Cb = - y j (w^)plog 2 (l + 7 )dw'd 7 , (27) 

0 0 


where p 


ve 


PsXl 


2dY 


iKijv) 


PsXi 


+ 


vKoiy) 


27(7+l-y^7(7+l)) 


and V = 2 


>Ci>C2-Ps Al A 2 


Proof: A proof is provided in Appendix D. ■ 

Although the TS factor does not appear directly in the ex¬ 
pressions in Proposition 3 and Proposition 4, the instantaneous 
CSI-based ctginr in (1251) has been employed in the derivations 
of ^out and C^. Therefore, the average throughputs in (I15bl) 
and (I15cl) for the SINR relay will be rewritten as 




Rr^ 


f E {(1 - asinr)R}, 7sinr > 7th 

\ 0? 7sinr ^ 7th 

^1(1 Ctginr) log2(l 3“ 7sinr)}'* 


(28a) 

(28b) 


/^sinr — 


2 sr^ 


"^SR CsR+^)9piD •Air I ho 7 


(23) 


which also satisfies the non-oscillatory condition in (O. Sub¬ 
stituting (1231) into ([T^ . we can express the e-SINR as 


7sinr 


_'^sr'^rd_ 

+2\/ 7sr ('Tsr+1)7rd •Air I ho 7 


(24) 


Obviously, fully utilizing the relay-harvested energy is not 
prerequisite in designing the relay gain /3sinr- Eor example, the 
redundant energy can be harvested by the relay in addition to 
the necessary energy harvested to support the relay gain /3sinr- 
To extend the relaying transmission time as long as possible, 
the TS factor is designed as small as possible such that the 


C. Target Relay 

When the SINR relay is employed, exact knowledge of 
the channel SNR 7 ^^,^ has to be exploited, which can be 
estimated only at the destination. A feedback channel is then 
required if Ofsinr is computed locally at the relay. Therefore, 
this subsection proposes a simplified relay control scheme, 
i.e., the target relay, that aims to achieve a target e-SINR 7. 
To avoid using knowledge of the target e-SINR should 
satisfy 7 < Denoting the e-SINR achieved by the target 

^The numerical computation of this expression can be quickly performed 
in softwares such as Matlab and Mathematica. 
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relay as 7tar (7tar = 7 ), the target relay is designed to 
maximize its e-SINR, i.e., to design the relay gain /3tar such 
that 7tar = 7sinr. Denoting the TS factor for the target relay 
by atar and substituting /itar — i^^to ([TtI) and (fTSl) . we 

can write the relay gain and e-SINR as 


and 


Ptar 


_ A^tar7sR _ 

^+^SR 7 sr 1 ^0 P 


(29) 


= _ /^tar7sR7RD _ /QQX 

7sR+(Mtar7sR“^r-|/ioP + l)(Attar7RD+l) ’ 

respectively. For a given 7, by eliminating from the 
equation pair {7tar = 7,7sinr = 7 }, the TS factor is given 

by (7tar = 

_(7sr + ^)(7sr~7)_ 

(7sR + l)(7sR~7+''777sR*^’"l^oP)-|-'^7sR‘^’"l^oP'y/7(7+l)7sR (7 sr + 1) 

Since o^tar harvests only a necessary energy to support the 
relay gain f3tar, the relaying transmission time (1 — atar)T 
is longer than those of other as satisfying a > atar- Also, 
since 7 < atar satisfies 0 < atar < 1- When 7 > 
we have o^tar < 0. In this case, no time is assigned for EH 
such that information relaying fails due to a lack of power. 
Alternatively, we can reset TS factor by atar = 1 such that 
only EH is implemented for the entire time block. 

Proposition 5: The outage probability achieved by the tar¬ 
get relay is given by 


^out — 1 ~ 


00 

a7 / 


l+K 




~2 — ^ 

ip 


1 + K-Ar 


(31) 


where p = 


_d7d^7tha>_ 

>^2A2((c^ -l)7th-l)((c^-l)>^lPs^-C^r^r) 


and uj = 


/ (7+l)(d7^o-^+£iPa2:) 

V CiPs^z 

Proof: A proof of this proposition is similar to the proof 
of Proposition 3. ■ 

For the target relay, the effective EH and relaying trans¬ 
mission fail when < 7 . In this case, the e-SINR 7 tar 
does not exist because of an impractical /itar- Although the 
derivative of Pout in Proposition 5 can still be obtained by the 
mathematical manipulation, it cannot be used to represent the 
PDF of 7 tar because of the discontinuity of 7 tar- Therefore, 
finding the PDF of 7 tar to evaluate = E{log 2 (l + 7 tar)} 
is difficult. As an alternative, the ergodic capacity achieved by 
the target relay will be investigated by simulations. Similar to 
the SINR relay, the average throughputs for the delay-limited 
and delay-tolerant transmissions are given by 


TO _ f 0 ^tar)P }5 Ttar ^ 7th 

"1 0, 7tar < 7th ’ 

^DT = E{(1 - atar)log 2 (l +7tar)}. 


(32a) 

(32b) 


D. Direct Transmission 

Due to an energy loss resulted from the EH circuit and 
the two-hop radio propagation, SWIFT FDR suffers from a 


degradation in the EE. However, since SWIFT FDR provides a 
relative freedom besides the source-destination direct transmis¬ 
sion, it results in an additional performance advantage. In this 
subsection, we consider the source-destination transmission 
scenario with a direct link, where an SWIFT FDR node can 
be employed to assist the information transfer. 

For the transmission scheme exploring both the direct link 
and SWIFT FDR, we assume that the signal transmitted from 
the source directly and the delayed replica transmitted from 
the relay are fully resolvable by the destination, so that they 
can be appropriately co-phased and merged via maximum ratio 
combing (MRC). The received e-SINR at the destination can 
be expressed as 

TmRC TsD 3“TsRD5 (33) 

where 7 g^ = is the SNR achieved by the direct 

transmission, hs is the channel coefficient of the direct link, 
and ds is the distance between the source and destination. 
Then, the instantaneous throughput of the direct transmissions 
with SWIFT FDR can be expressed as 

Rj =alog2(l + 7sD) + (l-a)log2(l+7MRc)- (34) 

It is easy to show that the outage probability of the direct 
transmission without SWIFT FDR is given by 

= l-e , ( 35 ) 

where A 3 is the mean of the exponential random variable \hs\‘^. 


Corollary 1: The outage probability of the direct transmis¬ 
sion with SWIFT FDR is upper-bounded by 

p(MRC) _ p(SD) p(SRD) 

^out,ub ~ ^out ^out 5 t^o; 

where the outage probabilities s achieved by the max¬ 

imum relay, SINR relay, and target relay are given by (l2Qab . 
(l26l) . and (|3T]) . respectively. 

Proof: A proof is provided in Appendix E. ■ 

Based on Corollary 1, the average throughput of the direct 
transmissions with the maximum relay in the delay-limited 
transmission mode is lower bounded by 

^^SO + SRO) = ^(1 _ p^SD) + (1 _ _ P77) )i?, (37) 

For the direct transmission with the SINR relay and target 
relay, where a depends on the instantaneous CSI, the average 
throughput can be evaluated by 


^SD+SRD _ E{ciQ(i? + C 2 (l - a)R}, 


(38) 


where ci = 


{J: 


Tsd ^ Tth 
7 sd < 7th 


and C 2 = 


1 , 

0 , 


7mrc ^ 7th 
7 mrc ^ 7th 


In the delay tolerant transmission mode, the ergodic capacity 
of the direct transmission witjiout SWIET FDR can be ex- 

pressed as (o, j , where r(-,-) 

in the upper incomplete gamma function 1391 Eq. (8.350.2)]. 

















TABLE I 

Simulation Parameters 


No. 

Parameter 

Value 

1 

Carrier frequency 

915 MHz 

2 

Bandwidth 

200 kHz 

3 

Fixed transmission rate R 

2 bps/Hz 

4 

Path loss at the reference distance 1 m: jC 

-30 dB 

5 

Distance between the nodes: di and d 2 

10 m 

6 

Path loss exponent m 

3 

7 

Means of dual-hop channel gains: Ai and A 2 

1 

8 

Noise power: and 

—95 dBm 

9 

Source transmit antenna gain 

18 dBi 

10 

Relay transmit antenna gain 

8 dBi 

11 

Rician iV-factor of the RSI channel 

6 dB 

12 

Energy harvesting efficiency 77 

0.8 


The throughput of the direct transmission with SWIET FDR 
can be expressed as 

^^SD+SRD) = + E{(1 - a) log 2 (l + 7mrc)}- (39) 

When only the statistical CSI is available, ([^ can be rewritten 
as + (1 - where ^ 

log 2 (l +7mrc) is ergodic capacity achieved by the MRC 
at the destination. Due to mathematical intractability, the 
throughputs expressed in (1^ and (l39l) and the related EE will 
be numerically investigated in the next section. 

E. CSI-related Issue 

Based on the pilot symbols sent from the source over dual¬ 
hop channels, CSI can be estimated to facilitate the wireless 
information and energy transfer (O, (TOl, ifTTIl . f\M , 1^ . 
Similarly to the works of Bdl and 1421 . a request-to-send 
(RTS)/clear-to-send (CTS) based channel estimation scheme 
can be employed at the beginning of the entire transmission 
and instantaneous CSI can be obtained. Surely, RTS/CTS 
based channel estimation incurs extra overhead and energy 
consumption. 

Moreover, the CSI estimation error may degrade the per¬ 
formance of SWIET FDR system. For the maximum relay, 
the bisection method needs full CSI, i.e., \hi\‘^ for i = 0,1, 
and 2 to evaluate the throughput. To compute a, the SINR 
relay also needs full CSI, while the target relay only needs 
I/lop and I/lip. The CSI estimation error can be modeled as 
hi = hi + hi for i = 0,1, 2, and 3, respectively, where the 
estimation error hi^e is a complex Gaussian random variable 
independent of hi. In addition, we assume that hi^e has the 
zero mean and the variance n\hi\^, where the scaling factor n 
is a positive real number representing the relative ratio between 
the estimation error and true CSI. In simulations, we apply the 
above procedure to investigate the effect of CSI estimation 
error on the performance metrics. 

IV. Numerical Results 

This section presents some numerical results to validate 
the analytical expressions developed in the previous section 


and discuss the EE performances for the considered relay 
control schemes. Unless otherwise stated, the parameters used 
in simulations are given in Table 1. 



(a) e-SINR versus Pg. 



(b) a versus Pg 

Fig. 2. Performance metrics in instantaneous transmission. 

Fig. 2 investigates the e-SINR and corresponding a versus 
the source transmission power in the instantaneous trans¬ 
mission. In Fig. 2, we focus on a single frame with the 
following channel setting: |/iop = 0.342, |/iip = 1.898, 
and |/i 2 p = 0.986. When the source transmission power Pg 
increases from 0 dBm to 50 dBm, the target e-SINRs for the 
target relay are set to increase linearly from 0 dB to 20 dB 
with a step size of 0.8 dB. As observed, the three relay control 
schemes achieve almost the same e-SINR and the highest e- 
SINR is achieved by the SINR relay. In the region of the 
low transmission power, the e-SINR achieved by the target 
relay is slightly lower than that of the maximum relay. When 
Pg increases, the e-SINR gap between the maximum relay 
and target relay becomes negligible. Fig. 2(b) shows that TS 
factor decreases for the three relay control schemes when Pg 
increases. This result shows that the relay node can harvest 
an enough energy in a short time when Pg is large enough 
{Pg > 20 dBm). In addition, since the SINR relay applies the 
largest a to harvest energy, this leads to the shortest relaying 
transmission time. From the point of view of the EE, Fig. 2 
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shows that the SINR relay achieves its EE depending on the 
e-SINR more than the effective relaying transmission time, 
whereas the target relay achieves its EE with its target e-SINR 
and corresponding (1 — atar)^- Note that the e-SINR and TS 
factor of the maximum relay are obtained by the bisection 
searching, which will result in the highest EE in instantaneous 
transmission, as depicted in Eig. 4. 


(a) 



Fig. 3. The instantaneous throughput under CSI estimation error. 

The impacts of CSI estimation error on the instantaneous 
throughput are illustrated in Eig. 3, where the channel re¬ 
alization is as same as that of Eig. 2. In Eig. 3(a), we set 
Ps = 30 dBm, 7 = 12 dB, and assume that ho, hi, and h 2 
suffer the same level of CSI estimation error. The curves in 
Eig. 3 are obtained by averaging over 20,000 random CSI 
estimation errors. As observed in Eig. 3(a), the throughputs of 
all the relay control schemes decrease when k increases. When 
K = 0 . 1 , the throughputs of the maximum relay and target 
relay become about 0.20 bps/Hz lower than that of perfect CSI 
(hz = 0), whereas the throughput of the SINR relay decreases 
about 0.27 bps/Hz from that of the perfect CSI. This result 
indicates that the SINR relay is more likely to be affected by 
the CSI estimation error than the maximum relay. Eor the SINR 
relay, Eig. 3(b) plots the throughput versus Ps under different 
KS. As observed in Eig. 3(b), the throughput of the SINR relay 
decreases slightly for /i: = 0.1 and hc = 0.2. When Pg increases, 
the throughput degradation also increases. Moreover, Eig. 3(b) 
also depicts the throughput degrading for k = 0 . 8 , which is 
corresponding to the scenario with a serious CSI estimation 
error. When Pg = 30 dBm, the throughput degradation is about 
1 bps/Hz for hi = 0.8. This result shows that the SINR relay 
can not work well in the case of a serious CSI estimation error. 

Eig. 4 illustrates the impacts of \ho\‘^ and Pg on the EE of 
instantaneous transmission. In Eig. 4(a), we set Pg = 30 dBm, 
7 = 12 dB, hi and h 2 are the same as those of Eig. 2. As 
a result, the EE achieved by the SINR relay is very close to 
that of the maximum relay. When |/ioP increases, the EEs of 
the maximum relay and SINR relay decreases simultaneously. 
Eor the target relay with 7 = 12 dB, its EE approaches to 
that of the maximum relay in the region of |/ioP between 
0.2 and 0.8. Beyond this region, the EE of the target relay 


X 10 ® 




0 10 20 30 40 50 

Ps (dBm) 

(b) EE versus Pg 

Fig. 4. EE in instantaneous transmission. 

degrades slightly compared to that of the maximum relay. The 
curves of EE versus Pg are plotted in Eig. 4(b), where we fixed 
I/lop = 0.1. Eor the target relay, 7 increase from 0 dB to 20 
dB with a step size of 0.8 dB. In the considered whole region 
of Ps , the EEs of all the three relay control schemes are at the 
same level. In addition, the EEs of all the three relay control 
schemes decrease when Pg increases. Actually, the maximum 
EEs are achieved by the three relay control schemes when Pg 
is less than —8 dBm for this channel realization. Since the EE 
maximization is a concave fractional program of Pg, the EE 
becomes a monotonically decreasing function of the considered 
practical Pg (Pg > 0 dBm). The results in Eig. 4 also show 
that the target relay is competitive compared to the maximum 
relay and SINR relay since the target relay has the closed-form 
TS factor and does not need the second-hop CSI. 

Eig. 5 shows the impact of 7 on the EE for the target relay 
under the same channel realization of Eig. 4(b). As observed, 
the maximum throughput can always be achieved by the target 
relay. Eor example, when Pg = 30 dBm, the channel SNR of 
the first-hop link is 7 sr = 85.78 dB, whereas the target e- 
SINR that achieves the maximum EE is about 7 = 15 dB, 
which satisfies 7 < 7 sr (Note that we set 7 = 12 dB when 
Pg = 30 dBm in Eig. 4(b)). Moreover, an approximately 10^ 
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bits/Joule EE decreases when 7 moves by approximately 4 dB 
away from 7 = 15 dB. This EE degradation is about 10% 
(within the same order of magnitude) to the EE achieved by 
the maximum relay. When the second-hop CSI is not available, 
the target relay is suitable for scenarios in which the priority 
is EE maximization. 



Fig. 5. EE versus 7 . 

The outage probability versus the transmission power and 
average RSI channel gain is depicted in Eig. 6 . In Eig. 6 (a), 
we set (Jo = 0.1 and 7 increases from 0 dB to 15 dB with 
a step size of 0.6 dB. Eig. 6 (a) shows that the analytical and 
simulation results match well. As observed in Eig. 6 (a), the 
outage probability gap between the maximum relay and SINR 
relay becomes larger when Pg increases. The reason for this 
phenomenon is that the time of energy harvesting of the SINR 
relay is a little more longer than that of the maximum relay, 
so that the corresponding e-SINR is more larger than that of 
the maximum relay in the region of high Pg. Note that the TS 
factor of the SINR relay is just a little larger than that of the 
maximum relay, which will not degrade the EE dramatically, as 
verified by the results of Eig. 7. In Eig. 6 (b), we set Pg = 30 
dBm and 7 = 8 dB. As observed in Eig. 6 (b), the outage 
probability of the maximum relay varies very slowly when erg 
is smaller than 0.1. However, when the value of Cq is smaller 
than 0.1, the outage probability of the SINR relay and target 
relay decreases dramatically with a decreasing Moreover, 
Eig. 6 (b) also verifies that the SINR relay achieves the smallest 
outage probability. 

Eig. 7 shows the EEs of all the three relay control schemes 
for delay-limited transmission, where Pg, and 7 are the 
same as those of Eig. 6 . As a result, the SINR relay achieves 
a highest EE in the whole region of Pg. This result verifies 
that a larger e-SINR and a proper aginr of the SINR relay 
outperforms the maximum relay in both the throughput and 
EE. Eig. 7(a) also shows that the EE gap between the SINR 
relay and maximum relay in the low region of Pg is larger than 
that in the high region of Pg . In addition, in the low region of 
Pg, the EEs of all the relay control schemes increase with the 
increasing of Pg, whereas when Pg reaches a threshold, the 
EE begins to decrease with the increasing of Pg. As observed 
in Eig. 7(a), the larger EEs are obtained in the low region 



(a) Outage probability versus Ps 



(b) Outage probability versus (Tq. 


Fig. 6 . Outage probability in delay-limited transmission. 

of Pg, e.g., from 5 dBm to 15 dBm. Thus, a relatively high 
outage probability achieved in the low region of Pg becomes 
preferable from the point view of EE. Eor example, Eig. 6 (a) 
shows that the SINR relay achieves Pout = 0.8 when Pg = 10 
dB, whereas Eig. 7(a) shows that the SINR relay achieves the 
highest EE when Pg = 10 dB. Although the throughput is a 
monotonically non-decreasing function of Pg (not shown in 
Eig. 7), the results of Eig. 7(a) reveal a different changing 
trend for the EE versus Pg. Eig. 7(b) investigates the EE versus 
(Tq with a middle transmission power, i.e., Pg = 30 dBm. As 
depicted in Eig. 7(b), the SINR relay achieves the highest EE in 
a large region of a^. Only when < 0.03, the EE of the SINR 
relay is just a little lower than that of the maximum relay. Eig. 
7 also shows that the target relay achieves a competitive EE 
performance. Although the EE of the target relay is lower than 
that of the maximum relay, the EE gap between the maximum 
relay and target relay decreases with the increasing of Pg. 

The ergodic capacities are examined in Eig. 8 , where erg = 
0.4, Pg and 7 are set as same as those of Eig. 7. Eor the ergodic 
capacities of the maximum relay and SINR relay, the analytical 
expressions match well with the simulation results. The ergodic 
capacity achieved by the SINR relay is the highest among all 
the three relay control schems since the SINR relay always 
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(a) EE versus Pg. 



(b) EE versus ctq 


Fig. 7. EE in delay-limited transmission. 

achieves the highest e-SINR. Note that the simulated ergodic 
capacity of the maximum relay is achieved with the bisection- 
obtained a and only statistics of CSI have been employed 
in the bisection searching. Fig. 8 also shows that the target 
relay achieves a competitive ergodic capacity compared with 
the maximum relay. 

Fig. 9 examines the EEs for the delay-tolerant transmission. 
In Eig. 9, Ps, (Tq, and 7 are set as same as those of Eig. 8 . Eig. 
9(a) shows that the EE gap between the maximum relay and 
SINR relay is very small. In the middle and high regions of Pg 
(Pg >4.5 dBm), the SINR relay achieves an EE higher than 
that of the maximum relay. Moreover, in the low region of Pg 
(Pg <4.5 dBm), the SINR relay achieves an EE lower than 
that of the maximum relay. Eig. 9(a) also shows that the EEs 
of all the three relay control schemes monotonically increase 
with the decreasing of Pg for the considered region of Pg . Note 
that the maximum EEs of all the three relay control schemes 
are obtained by further decreasing Pg into the extremely small 
region (around —20 dBm), so that such a small transmission 
power is out of practical use. Eig. 9(b) shows that the EE of the 
maximum relay and SINR relay decrease with the increasing 
of (Tq. Eor the considered target e-SINR, Eig. 9 also shows 
that the target relay achieves a competitive EE compared with 


the maximum relay. 



(a) Ergodic capacity versus 



(b) Ergodic capacity versus (Tq. 

Fig. 8. Ergodic capacity in delay-tolerant transmission. 

Eig. 10 investigates the impact of the relay position on the 
system performance. In Eig. 10, we set Pg = 30 dBm, = 
0.1, ds = 20 m, and di P d 2 = dg. As observed, the highest 
outage probabilities are achieved by all the three relay control 
schemes when di/d^ = 0.5. When the relay moves away 
from the middle position, the outage probabilities decrease. 
Eor the delay-tolerant transmission, the lowest throughputs are 
obtained when the relay is placed midway. The reason for 
this is that the e-SINR satisfies oc {did 2 )~^ and 73 ^^^^ 
achieves its smallest value when di = ^ 2 . Moreover, Eig. 10 
shows that the lowest EE is obtained when the relay is placed 
midway. Therefore, the relay should be placed closer to the 
source or destination to achieve a better system performance. 

Eig. 11 compares the system performances of the direct 
transmission with and without SWIET EDR in the delay- 
limited transmission mode. In Eig. 11, we set erg = 0.1, 
di = d 2 = 10 m, ds = 20 m, and consider R = 2 bps/Hz 
and 6 bps/Hz, respectively. Since the EEs of all the three 
relay control schemes are at the same order of magnitude, 
we consider only the SINR relay to clearly show the curves. 
As observed in Eig. 11(a), the outage probability becomes very 
small when the direct link is available. In the low region of Pg, 
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(a) EE versus Pg. 



(b) EE versus ctq 

Fig. 9. EE in delay-tolerant transmission. 

the outage probability of the direct transmission with SWIET 
FDR is almost the same as that of direct transmission. After 
Pg increases to a threshold value, the outage probability of the 
direct transmission with SWIET FDR can be further decreased 
compared to that of the direct transmission. However, when R 
increases from 2 bps/Hz to 6 bps/Hz, the outage probability 
decreasing requires a higher Pg, which results in a great 
EE decreasing, as depicted in Fig. 11(b). Since the outage 
probability of the direct transmission is already very small. 
Fig. 11(a) suggests that SWIET FDR is benefit in assisting the 
low-rate and high quality-of-service transmission. For example, 
when Pg = 20 dBm, the outage probability of the direct 
transmission is about 1.2 x 10“^ for R = 2 bps/Hz, which 
can be further decreased to 4.0 x 10“^ with the assistance 
of the SWIET FDR. Another important observation is that 
the EE of the direct transmission (with and without SWIET 
FDR) is monotonically decreasing when Pg increases. As a 
result, the EE gap between the direct transmission and relaying 
transmission decreases when Pg increases. For R = 2 bps/Hz, 
the EE gap almost disappears when Pg surpasses 30 dBm. 
The reason for this phenomenon is that the contributions of 
different small outage probabilities (e.g.. Pout < 0.1) to the 
EE become relatively small when these outage probabilities 


are obtained by the relatively large PgS. Therefore, it can 
be concluded that the SWIET FDR is an energy-effective 
scheme for the low-rate transmission when the direct link is 
not available. 


(a) (b) 



di/ds 

Fig. 10. Impact of relay position on system performance. 

Fig. 12 compares the system performances of the direct 
transmission with and without SWIET FDR in the delay- 
tolerant transmission mode. In Fig. 12, the simulation pa¬ 
rameters are the same as those of Fig. 11. As observed, the 
ergodic capacity of the direct transmission is much higher 
than that of relaying transmission. When Pg = 30 dBm, 
the ergodic capacity gap between the direct transmission and 
relaying transmission is about 17.7 bps/Hz. More over, when 
SWIET FDR is employed to assist the direct transmission, 
the ergodic capacity almost keeps the same as that of the 
direct transmission, since the e-SINR enhancement due to 
SWIET FDR is relatively small compared to that of direct 
transmission. As a result, the EE of the relaying system 
assisted direct transmission is almost the same as that of the 
direct transmission. Thus, the performance enhancement for 
the delay-tolerant direct transmission by employing SWIET 
FDR is trivial. 

V. Conclusion 

This paper has studied three relay control schemes to 
improve EE of SWIET FDR systems. The EE maximization 
problem has been formulated as a concave fractional pro¬ 
gram of source transmission power, so that the three relay 
control schemes are designed separately. Under Rician fading 
RSI channel, analytical expressions of outage probability and 
ergodic capacity have been derived for the maximum relay, 
SINR relay, and target relay. For the maximum relay, the EE 
maximization problem is shown to be a concave function with 
respect to TS factor, so that the optimized TS factor can be 
obtained by the bisection method. By employing instantaneous 
CSI, the SINR relay and target relay with collateral TS factors 
have gained advantages in improving EE. In delay-limited and 
delay-tolerant transmissions, numerical results show that the 
SINR relay achieves a better performance than that of the max¬ 
imum relay in terms of outage probability, ergodic capacity. 
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and EE, whereas the target relay achieves a competitive system 
performance without requiring CSI of the second-hop. It has 
shown that a worse outage performance may result in a higher 
EE in delay-limited transmission. In low-rate delay-limited 
transmissions, numerical results also show that SWIET EDR 
can achieve a competitive EE compared to source-destination 
direct transmission. When the relay is placed midway between 
the source and relay, the worst EE is achieved, as well as the 
outage probability and erogdic capacity. 



(a) 



Fig. 12. Performance in delay-tolerant direct transmission. 


number. Pout can be simplified as 


Prf/i 2 p< 



(i\hi p-t-6 

c\hi\^-d\hi\ 

1 , 

1 , 


l^ol^ < ^ and < f. 
1 ^ 01^10 


(a) Outage probability versus Ps 



(b) EE versus Pg. 


Eig. 11. Performance in delay-limited direct transmission. 


Appendix A: A proof of Proposition 1 


Substituting ([Ull into Pout = Pr(7max < 7th), the outage 
probability is given by 


out — 


= Pr ( 




= Pr \h 


(I 


7sr+(/^ 7SR^r I ^0 P + 1) (iU7RD+1) 
2 ^ q|^iP+^ 


<7th 


< ^ 




)■ 


(A.1) 


Denote the PDF of |hip and the cumulative distribution 
function (CDF) |/i 2 p by and = 

Pr(|/i 2 p < z) = 1 — e respectively, where Ai and A 2 are 
the means of |/iip and |/i 2 p, respectively. Then, Pout can be 
evaluated as 


00 00 

Pout= j j 

—-— d 
f^'yth c 

+ / / 4o|7+/|.i|2+d^du; 
w=0 z=0 
1 

f^'yth 00 

+ j j /|.0|7++i|2(7Pr(|/i2|" < ^f^)d^du; 

W=0 z=- 

c 

1 

A^^th 00 

= 1 — ^ J J f\h |2('^)^ ( + (cz‘^-dz)X2 ) dzdw, (A.2) 

w=0 z=- 


where a = CiPgd'^d'^al^thil + pAr\ho\‘^), b = 
c = C\C 2 PsP{l - /iTthAl^ioP), and d = 
ClC 2 Psd'^( t‘^ fJL^ith- Given that the term c\hi— d\hi\‘^ can be 
positive or negative and | /12 P is always greater than a negative 


where a = CiPsd'^d'^a‘^^th{l + pArw), c = C\C 2 PsP{l — 
PlthArw), and is the PDE of |/ioP given by ([T9l) . 

The analytical expression for Pout in (IA.2I) cannot be further 
simplified. However, at the high SNRs, has an approx- 
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imation as —so that Pout can be approximated as 


out^ 


1 - 




w=0 

c 

1 _ 

A^T'th OO 


A 1 

a:;=c2— 


-— [ If, ,o(w)e (cAi+a;A 2 )g cAi 

cAi J J 


dxdri; 


w=0x=0 


where the infinite-series representation of Io{z) (391 Eq. 
(8.447.1)] has been applied in the step (a). Now, the first term 
in the RHS of (IB.lb can be shown similarly as 

E{log2(l + M7rd + 


CXJ 

ThJ 4„|2Hg1;2 ( 


yCl >^2 Ai-Ps A1 A 2 (l + ^iwArW) OjOjldA 

d^d^al 1,0 ) 


dw. 


Then, the desired result follows immediately. 


j f^,,^^2iu))pKi{p)e -Adw, 


(A.3) 


Ce = E{log2(l + M7 rd + Ai^7RD A|/io|^)} 


-E{log2(l + 

p|/i 2 p, where the 
^0 Given \ho\^^ 

side (RHS) of (ED can b( 

E{log 2 (l+ Ai^ 7 RDA|(io|^)| \hof} 




Appendix C: A proof of Proposition 3 

Substituting (l24b into Pout = Pi*(7sinr < 7th), the outage 
probability is given by Pout = 


where p = y 0 X^X 2 ^i(') first-order modified 

Bessel function of the second kind (391 Eq. (8.432)]. The last 
equality in (IA.3b is obtained by applying -^^dx = 

EH Eq. (3.324.1)]. 

Appendix B: A proof of Proposition 2 

In the RSI dominated scenario, the ergodic capacity achieved 
by the maximum relay can be expressed as 


\TAj\h |2 / a|hi|^+6+CY^|hi|2(£iPd^iP+c^) \ |2 ^ c^rTthC^? 

<-z;(ZTft^i7F:w?- } 1 

1 , 

where d = CiPsAr\ho\'^dfd^ajjth{'i- + ‘^Jth), 
b = c = 2Ar\ho\'^d'^df(7l'yth 

^CiPsil + 7 th) 7 th,^ and d = df-a^. Substituting 
/ 2 (- 2 ^) — / 2 (- 2 ^) of (fTgl) . and F (z) = 

Pr(|/i 2 p < 2 ;) = 1 — into the above expression of Pout, 
the outage probability can be written as 


(B.l) 


Define x = |/ii P |/i 2 P, where the PDF of x is given by f{x) = 
Given |/ioP, the second term in the right 
hand side (RHS) of (IB.lb can be evaluated as 


Pout — 1 ~ X 


^ ^ w (^az-\-h-\-c^y z{CiPsZ-\-d)^ 

//|^q| 2 ('^)c i: 2 A 2 (/:iPs^- 7 thC«)^ dwdz 


JO^lPs 




_2 KAro-Qp 

— 1 L- I _ l-\-K _ Ai l + K + Arcx'^P tC D 

Xi I l+K+Ar-a^p^ 


^iPs 

where a = CiPsd^dfahthAA + 27 th), b = 
7j\A, c = 2df d^(T^ 7 thA\/EiPs(l + 7th)7th, and 


(a) 


Gli (a)go(2y7g) dx 

0 P 


(b) 


1X2^4,2 


£i£2Af^P3AiA2>4^|hop 

dTdX^a'j, 


1,0 ) 


(B.2) 


>^2 A 2 (>c 1 Ps 2 ; - 


-. (C.2) 


where G^^^(x) is the Meijer G-function (39l Eq. (9.301)]. In 
( 1121 ), we have used the relationship (431 Eq. ( 8 .4. 6 .5)] in the 
step (a) and the integral identity (39l Eq. (7.821.3)] in the step 
(b), respectively. Then, the second term in the RHS of (IB. lb 
can be expressed as 

E{log 2 (l + 77RD4r|/lo7} 

OO 

1 /" r ( C1C2APs^i^i-Art o,o,l,l^ 

1112 J /|,„|2(OW,2 dTdjvl 1,0 jdf 


(a) e-K 


if In 2 


E 

n=0 


dY^d'^c 7 ^K(K + l) 

ClC2P‘^ Ps^l>^2-^rPo 

(«!)" 


This proves Proposition 3. 

Appendix D: A proof of Proposition 4 
The ergodic capacity can be evaluated by = 

/o“ / 7 =inr( 7 )log 2 (l + 7)d7, where is the PDF of 

7 sinr- By replacing 7 th with 7 in (1C.lb . the CDF of 7 sinr 
can be obtained. Then, the PDF of 7 sinr can be computed by 
fisinril) — However, the complicated expression of 

p in (10211 leads to a huge expression for such that 

the numerical evaluation of becomes difficult. Fortunately, 
a high SNR approximation can be applied to simplify the 
expression. By using the similar procedure as in (IA.3b . the 
PDF of 7 sinr can be computed as 


^ pOC 

X(7^4 f d7d^cr^(K+l) I -1-n ^ /o 3\ / (7) = — / f\ui2(w)pdw, 

>^^ 1,4 £l£2Af2p^AlA2^rCr2 I 0 ,-l-n,-l-n,-ny 5 Jls.nrKIJ J\no\ V 


(D.l) 
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2d^ cr^7 

where p ^ 


and V is given by 


V 


^ ^ I ArwdT'd^alj(2j+l+2^/r{r-^^ 

~ ztz^piata; ■ 


Then, the ergodic capacity is given by 

^oo />oo 


poo poo 

Ce = - / / f\ho\^{w)plog2{l + 'j)dwd'j, 

Jo Jo 


(D.2) 


where the term p can be further evaluated as 


2dfa^vKi(v) 

P= - ^ I . ^ 


P.Ai 


+ e 


■ d[vKi{y)] 
d'y 


2d^al^ 


ve 




PsXl 


+ 


vKq{v) 


27(7+l-y^7(7+l)) 

Appendix E: A proof of Corollary 1 


(D.3) 


The outage probability of multiple relays assisted transmis¬ 
sion can be written as 


p(SD + SRD) 

-CduI 




(Tsd ^ Tth Tsrd) 5 Tsrd Tth 


0 , 


Tsrd ^ Tth 


(E.l) 
P by 


Denoting the CDF and complementary CDF of |/i 3 

and respectively, and denoting the PDF 

7srd by tbe outage probability in (IE.II) can be 

expressed as 


p(SD + SRD) 

-CduI 



/•7th 


/ F J 


Jo 

= 

p(SRD) _ f 
-^out / 


Jo 

(a) 

p(SRD) f _ 
-^out 

< 

p(SRD) p 
out ^1^3,2 

= 

p(SD) p(SRD) 
-^out -^out 

) _ 



(7th - 7) A 

'Tth 


(7)d7 


^|,3|2(7th-7)/7 


(7)d7 


\hs\‘^ 

(7th) 


(7th - 7c 


(E.2) 


achieved by the SWIFT FDR and the constant jc C [0,7th]. 
In the step (a) of (IE.2b . we have applied the weighted mean 
value theorem ||44l Theorem 3.16]. 
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